Unquenched calculation of as from Green functions: progress report.* 
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We present preliminary results on the computation of the QCD running coupling constant in the MOM 
scheme and Landau gauge with two flavours of dynamical Wilson quarks. Gluon momenta range up to about 
7 GeV (/3 = 5.6, 5.8 and 6.0) with a constant dynamical quark mass. This range already allows to exhibit some 
evidence for a sizeable correction to the asymptotic behaviour, as in the quenched approximation. We 

find A^"^ = 264(27)MeV x [a"^ (5.6, 0.1560)/2.19 GeV], which leads to as{Mz) = 0.113(3)(4). In view of the 
systematics error to be controUed, this encouraging result is more a preliminary indication than a real prediction. 



1. Introduction 

From the systematic study of as in the 
pure Yang-Mills case |jl| by the Green func- 
tions method we retain that the "asymme- 
tric" 3-gluon vertex (i.e. with incoming momenta 
(/x^,0,/L<^)) is a reliable tool. We also recall that 
power corrections were exhibited and not negli- 
gible up to 10 GeV they were eventually 
traced back to an {A^) condensate through an 
OPE study 1). 

This whole work is now undertaken with two 
flavours of dynamical quarks [^. Following the 
conclusions of and denoting by as, pert the 
perturbative running coupling up to four loops, 
we fit the data according to 



MOM 



(1) 



on a large energy window obtained by combining 
several lattice settings corresponding to the same 
renormalised dynamical- quark mass expressed in 
physical units. Starting from a calibrating lat- 
tice setting (/?, Ksca)=(5. 6, 0.1560) taken from [||, 
we search Ksca's to fulfill this mass condition for 
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and 6.0. To this aim, we perform an explo- 
ration of the space of bare parameters (/3, Ksoa)- 
After that we are in position to fit the data ac- 
cording to eq. (0) and get A^^-^ At last, we 
derive a preliminary estimation of ajjsiMz) and 
briefly discuss uncertainties on this result. 

2. Exploration of parameters space (/?, Ksca) 

Given a calibrating lattice setting, we only need 
to evaluate lattice spacing ratios between two dif- 
ferent sets (/9, Ksea); this is done by ensuring the 
smoothness of as on the whole energy window 
under consideration. Quark masses are estimated 
using the axial Ward identity with Kgoa = '^vai 

We have taken Z^=0.77(l) and Zp=0.54(l) from 
101 in the RI-MOM scheme at 3 GeV, and con- 
verted niq to MS This method provides us 
with first estimates (see table H) which will be 
refined later on. We see that the masses cor- 
responding to the lattice settings (5.6,0.1560), 
(5.8,0.1525) and (6.0,0.1505) are almost equal. 

Because of the computational cost of this ex- 
plorative stage, we only worked on small volumes. 
However two checks indicate that the information 
we gathered is relevant to tune (/9, Ksoa) for runs 
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Table 1 

First estimates from 
computed in the MS 



our runs; 
scheme at 3 



the masses are 
GeV. 



/3 


^sca 


Volume 


a-' (GeV) 


m^ea (MeV) 


5.6 


0.1560 


24716^ 


2.19(8) 


164(7) 


5.6 


0.1575 


16^ 


2.42(9) 


79(3) 


5.8 


0.1500 


16'* 


2.45(13) 


325(18) 


5.8 


0.1525 


16-* 


2.76(7) 


173(4) 


5.8 


0.1535 


164 


2.91(18) 


103(16) 


5.8 


0.1540 


16" 


3.13(13) 


64(5) 


6.0 


0.1480 


16" 


3.62(10) 


391(12) 


6.0 


0.1490 


16" 


3.73(13) 


308(12) 


6.0 


0.1500 


16" 


3.78(14) 


213(3) 


6.0 


0.1505 


16" 


3.84(15) 


169(8) 


6.0 


0.1510 


16" 


3.96(16) 


96(4) 



on larger volumes. Firstly we did not observe any 
significant differences on as at large momenta 
between our two runs at (/3, Ksoa)=(5. 6, 0.1560). 
Secondly there could be a deconfinement and / or 
chiral phase transition at very small volumes, in- 
validating the relation mp oc (m^ + m^), where 
mp denotes the mass of the lightest pseudoscalar 
fitted on our lattices, or, strictly speaking, the 
pion mass mixed (due to the small lattice length) 
with some excited states. We found empirically 
that 



= 2Br71sca + y 



(3) 



with B = 2.74(5) GeV and r = 1.41(5) GeV^fm" 
from a best fit with a x^/d.o.f .=0.57. This re- 
lation indicates that we suffer from strong but 
smooth finite volume effects on mp (which is 
an IR-sensitive quantity). We have neverthe- 
less some control on it, since at the continuum 
limit at {f3, Ksoa)=(5. 6, 0.1560) we recover rn|, ^ ~ 

0.432(18) GeV^, in good agreement with SESAM 
Q; furthermore, we get (m„ -|- md)/2 ~ 3.6 MeV 
from the pion mass and rris ~ 90 MeV from the 
kaon mass. This compares fairly well to other 
lattice estimates. 



Fitting 

^ MS 



and power corrections 



Keeping the set of lattice settings with constant 
dynamical mass, we refine our estimations of lat- 



tice spacings by considering a polynomial fitting 
all the lattice data (see figure |l|) except few points 
corresponding to n = (L^)^/(47r^)<2 — 4 (where 
L is the length of the lattice), presumably affected 
by strong finite volume effects ||l|. The stability 
of the fit is fairly good and yields to 

0-1(5.8,0.1525) 
0-1(6.0,0.1505) 



= 2.85 ±.09 ±.04 

a~^(5. 6, 0.1560) 



^ 2.19 GoV 

3.92 ±.11 ±.07 

^(5.6,0.1560) 



GeV 



(4) 



X 



2.19 GoV 



GeV 



The central value corresponds to a cut at n > 3, 
the first error is statistical and the second is sys- 
tematic. We are going to use these sharper values 
from now on. Since the error on the calibrating 
lattice spacings propagates trivially, we drop it 
till the discussion of systematics. 




Figure 1. Matching of lattice spacings. 



Having fixed the lattice spacings, we fit A^-^ 
according to e(7.(|l|) on a wide energy window 
ranging from 2.6 GeV, giving (see figure 0) 



MS 



264(27)- 
2.7(1.2) 



-1(5.6,0.1560) 

2.19GeV 
a-i(5.6, 0.1560) 



2.19 GeV 



MeV 



GeV 



(5) 
(6) 



with a reasonable x^/d.o.f. ~ 0.6. 

At this point, we should stress that another 
fit without power corrections is also acceptable 
ixVd.o.f. = 0.96), giving A^=' = 345(6) MeV. 
Therefore more statistics is needed to determine 
the size of power corrections. 
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Figure 2. Asymptotic fit of as- 



4. Running of as to AIz and discussion 

Concerning finite spacing effects, let us first 
remark; that the use of non-improved dynamical 
quarks leads to 0{a) errors. This induces only a 
slight effect on the ratios of our small lattice spa- 
cings. The dominant error thus comes from the 
calibrating lattice spacing and propagates mul- 
tiplicatively. From [ll| and we estimate this 
systematic uncertainty to be ~ 50 MeV on A^~^ ; 
its effect on as{Mz) is an error of ±.004. The 
smooth junction of as from different lattices show 
that other overwhelming effects are absent. 

We neglect of course power corrections in the 
extrapolation to Mz- If we perform the same 
analysis (i.e. according to e(7.(|l])) with quenched 
data, we can extrapolate our Nf = 0, 2 results 
to Nf = 3 at, say, the charm threshold. A stan- 
dard running procedure |^ , taking the charm and 
beauty thresholds at 1.3 GeV and 4.3 GeV, gives 



as{Mz) 

as,pert{Mr) 



0.113(3)(4) 
0.283(18)(37) 



(7) 
(8) 



The second error comes from the error on the 
calibrating lattice spacing. This result is 2a below 
the world average experimental result as{Mz) = 
0.119(2) and the result of ALEPH H as(Mj) = 
0.334(22). It is compatible with the recent result 
of the QCDSF-UKQCD coll. [|l3|; older resuhs 
using NRQCD were closer to experiment [|4|. 

Being 2a below the experimental value is very 
encouraging. Work is now in progres to check in 



more details the effects of our large quark masses 
and small lattices. This will allow us to refine our 
preliminary result on as{Mz) in a near future. 
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